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investigated and the fractional Ostrogradski's formulation is obtained. The 
fractional path integral of both simple harmonic oscillator with an acceleration- 
squares part and a damped oscillator are analyzed. The classical results 
are obtained when fractional derivatives are replaced with the integer order 
derivatives. 

' I. INTRODUCTION 

a 

Fractional calculus deals with the generalization of differentiation and 
integration to non-integer orders. Fractional calculus has gain importance 
especially during the last three decades 1-5 . A large body of mathematical 
knowledge on fractional integrals and derivatives has been constructed. Frac- 
tional calculus, as a natural generalization of classical calculus, has played 
a significant role in engineering, science, and pure and applied mathematics 
in recent years. The fractional derivatives are the infinitesimal generators of 
a class of translation invariant convolution semigroups which appear univer- 
sally as attractors. 

Various applications of fractional calculus are based on replacing the time 
derivative in an evolution equation with a derivative of fractional order. The 
results of several recent researchers confirm that fractional derivatives seem 
to arise for important mathematical reasons 5-21 . 



1 On leave of absence from Institute of Space Sciences, P.O.BOX, MG-23, R 76900, 
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The fractional variational principles represents an important part of frac- 
tional calculus and it is deeply related to the fractional quantization pro- 
cedure There are several proposed methods to obtain the fractional Euler- 
Lagrange equations and the corresponding Hamiltonians However this issue 
is not yet complectly clarified and it requires more further detailed analysis. 

Quantization of systems with fractional derivatives is a novel area in the- 
ory of application of fractional differential and integral calculus. Schrddinger 
equation was considered with the first order time derivative modified to Ca- 
puto fractional ones in 22 . In this case the obtained Hamiltonian was found 
to be non-Hermitian and non-local in time. In addition, the obtained wave 
functions are not invariant under the time reversal. The quantization of 
fractional Klein-Gordon field and fractional electromagnetic potential in the 
Coulomb gauge and the temporal gauge were investigated very recently in 23 . 

Recently, the fractional variational principles and the fractional Euler- 
Lagrange were obtained 24-25 . 

Even more recently, the fractional constrained Lagrangian and Hamil- 
tonian were analyzed 26-27 . The notion of the fractional Hessian 27 was in- 
troduced and the Euler-Lagrange equations were obtained for a Lagrangian 
linear in velocities 26 . Besides, the Hamiltonian equations have been obtained 
for systems with linear velocities 28 . The classical fields with fractional deriva- 
tives were investigated by using the fractional Lagrangian formulation and 
the fractional Euler-Lagrange equations were obtained in 29 . 

Non-local theories have been investigated in several physical problems 30 . 
During the last decade the non-local theories were subjected to an intense 
debate 31-35 . A Hamilton formalism for non-local Lagrangians was developed 
in 34-35 , an equivalent singular first order Lagrangian was obtained and the 
corresponding Hamiltonian was pulled back on the phase space by using the 
corresponding constraints 34 . It was shown the space-time non- commutative 
field theories are acausal and the unitarity is lost 36-37 . The fractional La- 
grangians and Hamiltonians are typical examples of non-local theories. 

For these reasons the fractional quantization of field theory is an inter- 
esting issue to be investigated. 

In this paper we analyze the fractional Hamiltonian quantization of non- 
singular systems possessing higher order derivatives. 

The plan of the paper is as follows: 
In section 2 the 1+1 classical dimensional field theory analysis of non-local 
theories is briefly reviewed and the fractional generalization of Ostrograd- 
ski's formulation is presented. In section 3 the path integral quantization 
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of the simple harmonic oscillator with an acceleration-squares part is ana- 
lyzed. Section 4 is dedicated to the fractional path integral formulation of 
the damped oscillator. Finally, Section 5 is dedicated to our conclusions. 



II. FRACTIONAL FIELD THEORY 
A. Classical non-local theory 

Let us start with an ordinary local Lagrangian depending on a finite 
number of derivatives at a given time, namely 

L(q(t),q(t),-,Q (n) (t))- (1) 

The next step is to consider a Lagrangian depending on a piece of the 
trajectory q(t, A) for VA belonging to an interval [a, b] 

L non (t) = L(q(t + A)), (2) 

where a, b are real numbers. Therefore a non-local Lagrangian was intro- 
duced. In this case the action function corresponding to ([2]) is given by 



S(q) = J dtL non (t) (3) 
and the Euler-Lagrange equation corresponding of ([3]) are given by 



/ 



*3H 

Equations (BJ should be understood as a functional relation to be satisfied 
by physical trajectories, i. e., a Lagrangian constraint. These functional 
relations define a subspace Jr of physical trajectories Jr C J, in the space 
of all possible trajectories 32 ' 34 . The crucial point is that there is no dynamics 
except the displacement inside the trajectory, namely 

q (t)^q(t + X). (5) 
Let us introduce now the dynamical variable Q(t, A) as follows 

Q(t,X) = q(t + \). (6) 

If we consider a field Q(t, A) instead of a trajectory q(t), such that 
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Q(t,X) = Q'(t,X), (7) 

where Q = dQM) and Q'(t, A) = 9 ®q^ we obtain a field theory in one 
spatial and one time dimension, namely a 1 + 1 dimensional formulation of 
non-local Lagrangians 32 ' 34 . 

The coordinates and momenta are suppose to have the following forms 



Q(t, A) = £ e m {X)q {m \t),P(t,X) = E e m (A)p M (t), (8) 

where 

{q {n \t),P {m) (t)} = 5 n m (9) 

and 

\ m 

e m (A) = — , e m (A) = (-d x ) m 5(X). (10) 
ml 

Therefore, the Hamiltonian for 1 + 1 dimensional field becomes 

H(t, [Q, P}) = J dXP{t, X)Q'(t, A) - L(t, [Q]), (11) 

where P denotes the canonical momentum of Q. The phase space is T * J 
together with the fundamental Poisson brackets 

{Q(t,\),P(t,\')} = 6(\-\'). (12) 

The functional L(t, [Q]) is defined as follows 

L(t,[Q}) = J d\5(X)C(t,X). (13) 

By using ([TBI the primary constraint arises as given below 

<j>(t, A, [q, P]) = P(t, A) - J da X (X, -a)e(t; a, A) « 0. (14) 

Here e(t; a, A) and x(A, —a) have the following definition 
. dC(t,a) . e(A) — e(o) 

where e(X) is the sigma distribution. The Euler-Lagrange equation is 
guaranteed by itself 

<P~ip = J da£(t;a,X). (16) 



B. Fractional Ostrogradski's construction 

Higher -derivatives theories 38,39 appear naturally as corrections to general 
relativity and cosmic strings 40 . Unconstrained higher-order derivatives pos- 
sess specific features, namely they have more degree of freedom than lower- 
derivative theories and they lack a lower-energy bound. A method how to 
remove all these problems was presented in 41 . It was observed that the non- 
local formulation translates into infinite order Ostrogradski's formulation 3 ' 35 . 

In this section, we would like to derive both the Lagrangian and the 
Hamiltonian formalisms for non singular Lagrangians with fractional order 
derivatives starting from the Hamiltonian formalism of non local-theories 32 . 
Let us consider the following Lagrangian to start with 

L(q,t) = L(t,q a ™), (17) 
where the generalized coordinates are defined as 

q Qm = aD? m x(t), (18) 

where m is a natural number. 

To obtain the reduced phase space quantization, we start with the infinite 
dimensional phase space T * J(t) = {Q(t, A), P(t, A)}. 

The key issue is to find an appropriate generalization of (fTOI) for the 
fractional case. As it was pointed out in 32,34 the coordinates and the momenta 
are considered as a Taylor series. Therefore, the first step is to generalize the 
classical series to the fractional case. A natural extension is to use instead 
of factorial the Gamma function. In this way we introduce naturally the 
generalized functions 42 instead of e m (A) and e m (A) given by (fTOl) . 

As it is already known several fractional Taylor's series expansions were 
developed 3,43 , therefore we have to decide which one is appropriate for our 
generalization. Since we are dealing with fractional Riemann-Liouville deriva- 
tives we choose the generalization proposed in 44 , namely 

oo 

Q(t,A) = J2 e am (X)q (a -\t), 

m=—oo 

oo 

P(t,X) = £ e a -(X) P(am) (t) 1 (19) 

m=— oo 

where 
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e Qm (A) = (A A j°" e«"(A) = D?»*(A - A ), (20) 

and a m = m + a, with < a < l.Here Ao is a constant. The coefficients in 
([19]) are new canonical variables 

{q [am \p {a ,)} = <&'• (21) 
By using (I2TT) we obtain that 

oo 

Y: e a -(X)e am (X')=5(X-X'), (22) 



m=— oo 



and 

/+oo 
rfAe Qm (A)e Q ,(A) = r™,. (23) 
-oo m am 

Therefore, e am (X) and e Qm (A) form an orthonormal basis. 

We stress on the fact that (|22|) and fl23|) involve the generalized functions 
and the relations have the meaning in the sense of generalized functions 
approach 42,44 . 

The fractional Hamiltonian is now given by 

oo 

H= Yl p am q am+1 - L(q°,q am ). (24) 

m=— oo 

The momenta constraints become an infinite set of constraints 

oo ar 

^^W-E.Dr^^0. (25) 
The fractional Euler-Lagrange equations are as follows 

a (26) 

An interesting property of the fractional series proposed by Riemann and 
discussed by Hardy in 44 is that when a m becomes integers the usual form of 
Taylor series is obtained. Therefore one should notice that for integer values 
of a m we have 

*-W- E l- Jt ) evF*wr °< (27) 



1=0 



6 



which is the definition of Ostrogradski's momenta 38 . 

In this case the Euler-Lagrange equation for original fractional derivative 
Lagrangian 26-30 is given below 

t> D >mrr°- (28) 

Now, from this equation, for integer values of a m we obtain the Euler- 
Lagrange equation for higher derivative Lagrangian 32,34 ' 38 , namely 

n / dL(t) 

g {-jt) mm = °' (29) 

The constraints fl27|) and (129]) lead us to eliminate canonical pairs 

{q ai ,p ai }(i>n). 

In this case the infinite dimensional phase space is reduced to a finite 
dimensional one. The reduced space is coordinated by T * J n = {q ai ,p Ql } 
with I = 0,1, ...,n — 1. The Hamiltonian in the reduced space is given by 

n-1 

H = p am q am+1 - L(q°, q am ). (30) 

m=0 

One should notice that the canonical reduced phase space Hamiltonian 
( I3"U|) is obtained in terms of the reduce canonical phase space coordinates 
{q ai ,Pa t } with I = 0,1, ...,n — 1. In this case the path integral quantization 
of filed system is given by 

K = J J] dq a ™dp a -e l ^ dt ^o pam « am+1 - H ». (31) 

m=0 

We observe that when a are integers, we obtain the path integral for sys- 
tems with higher order Lagrangians 32,45-46 . 

III. FRACTIONAL PATH INTEGRAL QUANTIZATION OF 
A SIMPLE HARMONIC OSCILLATOR POSSESSING 
ACCELERATION-SQUARES PART 



The classical Lagrangian to start with is given by 
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L = i(l + eV)x 2 -^V-ie 2 x 2 . (32) 
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The fractional generalization of (I32p has the following form 

L' = \(l + eV)( t D£n(f)) a - i W V - ~e 2 [ t D«( t D>(*))] 2 . (33) 

The independent coordinates are x(t) and tD®x(t) respectively. Let us 
denote pf = p x and p% = P( t D a x(t)) " The fractional canonical momenta are 38 

Pl d t n:x(t) tUffii 3 t Df,(t) J,ft d t r>l« x (ty [6V 

By making use of fl33l) we obtain the forms of the fractional canonical mo- 
menta as given below 

pf = (1 + eV) t D£r(t) + e\r>l a x(t), (35) 
p a 2 =-e\Bl a x(t). (36) 
Taking into account f[3"5"j) the fractional canonical Hamiltonian becomes 

H = p a lt D a a x{t) + P5tV?x{t) ~ L, (37) 

and after taking into account (1331) . (j35p and (I36p the fractional Hamiltonian 
has the form 



1 (r> a ) 2 
H = r[2pf*D^(t) - + wV(*) - (1 + eV)( t D>(i)) 2 ] (38) 

By making use of (|38p the fractional path integral is given by 

# = /"^(tD^^))^^^^/^? 1 ^^^* ^^)-^}. (39) 



IV. FRACTIONAL PATH INTEGRAL QUANTIZATION OF 
DAMPED HARMONIC OSCILLATOR 

The Lagrangian for this system in Ostrogradski's notations 38 takes the 
form 9 

L=\™£ + il<£ /2 -V(q ), (40) 

where 

q a - = t D% n x,n = 0,1,2. (41) 
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Here ao — 0, a.\ — |, a 2 = 1 and go = x, gi = x, gi = tD| x, g2 = x. 
The expressions for canonical momenta are 

Po = hx(i/2) + imx {3/2 ), (42) 
Pi/2 = mx. (43) 

By using fj4*0l) the classical Euler-Lagrange equation of motion read as 9 

dV 

m'x + 72 + — — = 0. (44) 
ox 

The canonical reduced Hamiltonian has the following expression 

H = ^i + qi/2Po ~ z 2 q2 ^ + v{qo) - (45) 

Therefore, the path integral representation for the above system analyzed 
is given by 



K = Jdpexpi J I gipi/2 - + i^qy 2 ~ V(q ) J dt , ( 10) 



where d\i = dq dp dqi/ 2 dpi/ 2 . 

The path integral representation for (IIBI) is an integration over the canon- 
ical phase space coordinates (go,Po) and (gi/2,Pi/2)- Integrating over pi/ 2 and 
Po, we obtain 

K = J dq dqi/2 expi J Qmg^ - V(q ) + ^1/2) dt - ( 47 ) 
Equation (14"T|) can be put in a compact form as follows 

K = J dqoe'I (k^ii~y(io))dt ^ e * f(%? /2 )<ft ( 48 ) 

After performing an integration over qi/ 2 (JHJ) becomes 

K = C J dq e ^^ m ^- v{qo) ) dt dt, (49) 
where C represents a constant. 

V. SUMMARY 
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The interest in fractional quantization appears because it describes both 
conservative systems and non conservative systems as well. The fractional 
quantization of field theory is not an easy task, especially when the fractional 
Hamiltonian is involved. The fractional derivatives represent the generaliza- 
tion of the classical ones and therefore some of the classical properties are 
lost e.g., the fractional Leibniz rule, the chain rule become more compli- 
cated than the classical counterparts. On the other hand, the fractional 
calculus represents an emerging field and it describes better various phe- 
nomena from several area of science and engineering. The fractional path 
integral formulation deserves further investigations mainly because the frac- 
tional generalization of the classical case is not yet complectly understood. 
Namely, for a system possessing second class constraints it is difficult to find 
the corresponding fractional generalization. In addition, there is no frac- 
tional formulations of the classical secondary or tertiary constraints because 
the fractional Hamiltonian is not a constant of motion. 

In this paper we generalize to the fractional case the non-local theories 
in one space and one time dimensions via the infinite Ostrogradski's for- 
malism. The classical Taylor series involved in this problem are convergent 
because of the properties of the Dirac's delta function. Namely, the co- 
ordinates and the corresponding momenta are defined as Taylor series and 
the Ostrogradski's canonical pairs fulfill the classical Poisson's brackets com- 
mutation relations. The generalization to the fractional case of all above 
mentioned results is not straightforward because there exist many formula- 
tions for the fractional Taylor series. However, a powerful tool in fractional 
field theory is to work to the Riemann-Liouville derivatives because of their 
important property of integration by parts. Therefore, in this paper we focus 
on the fractional Taylor series involving the Riemann-Liouville derivatives. 
We assumed that the fractional Lagrangian density has a compact support in 
the x-directions. In this work we have obtained the path integral quantiza- 
tion for fractional generalization of a 1 + 1 dimensional non-local field theory. 
The path integral formulation for the simple harmonic oscillator with an 
acceleration-squares part as well as for the damped oscillator are obtained. 
It is worthwhile to mention that the general expression for the path integral 
leads to the path integral representation for systems with higher order La- 
grangians. 
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